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FormCalc is a matrix-element generator that turns FeynArts amplitudes up to one loop into a Fortran code for 
computing the squared matrix element. The generated code can be run with FormCalc's own driver programs or used 
with other 'frontends', e.g. Monte Carlos. Major new or enhanced features in Version 4.1 are: treatment of external 
fcrmions, phase-space integration, code-generation functions, extensions for the MSSM, the HadCalc frontend. 



1. Introduction 



FormCalc Q is a package for the calculation of Feynman amplitudes based on Mathematica and FORM Q. 
Amplitudes generated by FeynArts 0] are first simplified analytically and then converted to a self-contained Fortran 
code for the computation of the squared matrix element. Currently, diagrams up to one loop can be simplified, and 
kinematics are supplied for 1 — > 2, 2 — > 2, and 2-^3 processes. 

The present article describes the features added or enhanced in version 4.1. They fall into several categories: 



• Weyl-van der Waerden (WvdW) spinor formalism for external fermion lines, 

• Multidimensional phase-space integration with Cuba 1.2, 
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• FeynHiggs interface, 

• Non-minimal flavour violation (NMFV) for the MSSM, 



• Various useful scripts, 



• The HadCalc frontend, by M. Rauch, for hadronic calculations. 



X 

^ , 2. Weyl-van der Waerden spinor formalism 

_ _ i 

Amplitudes involving external fermions have the form M. = c iFii where the Fi are (products of) fermion 

chains. The textbook recipe is to compute probabilities, such as \M\ 2 = Yl7j=i c i c i F*Fj, and evaluate the F*Fj 
by standard trace techniques: \(u\ T \v)\ 2 — (u\ T \v) (v\ f |ti) = Tr(r \v) (v\ f \u) (u\). 

The problem with this approach is that instead of n of the -Fj one needs to compute n 2 of the F*Fj. Since 
essentially n ~ (number of vectors)!, this quickly becomes a limiting factor in problems involving many vectors, e.g. 
in multi-particle final states or polarization effects. 

The solution is of course to compute the amplitude A4 directly and this is done most conveniently in the WvdW 
formalism . The implementation of this technique in an automated program has been outlined in Q . 

The FermionChains option of CalcFeynAmp determines how fermion chains are returned: Weyl, the default, selects 
Weyl chains. Chiral and VA select Dirac chains in the chiral (w+/w_) and vector/axial- vector (I/75) decomposition, 
respectively. The Weyl chains do not need to be further evaluated with HelicityME, which applies the trace technique. 

The WvdW method has other advantages, too: Polarization does not 'cost' extra in terms of CPU time, that is, one 
gets the spin physics for free. Whereas with the trace technique the formulas become significantly more bloated when 
polarization is taken into account, in the WvdW formalism one actually needs to sum up the polarized amplitudes to 
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get the unpolarized result. There is also better numerical stability because components of fc M are arranged as 'large' 
and 'small' matrix entries, viz. 

V ki + ik 2 ko — h J 

Cancellations of the form ko — k = \fk 2 + m 2 — y/k 2 for m <C k can be avoided (e.g. by observing that ko — k = 
m 2 /(ko + k) for on-shell k) and hence mass effects are treated more accurately. 



3. Multidimensional integration with Cuba 

FormCalc uses the Cuba library [6j for numerical integration of multidimensional phase-spaces. Cuba provides 
four integration routines: Vegas, Suave, Divonne, and Cuhre. All four have a very similar invocation and can thus 
be interchanged easily, e.g. for comparison. The flexibility of a general-purpose method is particularly useful in the 
setting of automatically generated code. The following table gives an overview of the features of the Cuba routines; 
specific details on the implementation and the actual usage of the Cuba routines are provided in Q and shall not be 
repeated here. 



Routine 


Basic method 


Type 


Variance reduction 


Vegas 


Sobol sample 


quasi-Monte Carlo 


importance sampling 




or Mersenne Twister sample 


pseudo-Monte Carlo 




Suave 


Sobol sample 


quasi-Monte Carlo 


globally adaptive subdivision 




or Mersenne Twister sample 


pseudo-Monte Carlo 


+ importance sampling 


Divonne 


Korobov sample 


lattice method 


stratified sampling, 




or Sobol sample 


quasi-Monte Carlo 


aided by methods from 




or Mersenne Twister sample 


pseudo-Monte Carlo 


numerical optimization 




or cubature rules 


deterministic 




Cuhre 


cubature rules 


deterministic 


globally adaptive subdivision 



Apart from some important bug-fixes, Cuba Version 1.2 includes the following new features: 

• Pseudo-random sampling using the Mersenne Twister generator has been added to all Monte Carlo algorithms. 
This is useful not only for comparison, but also because in high dimensions the numerator in the convergence 
rate for quasi-random samples, C(log d_1 n/n), becomes noticeable. 

• Vegas can memorize its internal grid for subsequent invocations, to speed up integration on similar integrands. 

• Vegas can save its internal state in a file such that the calculation can be resumed e.g. after a crash. 

• There exists a one-stop invocation for the Cuba routines, 

subroutine Cuba(ndim, integrand, result, error) 

with which the number of parameters that need to be passed is reduced to a minimum. (Note that this 
particularly simple form is special to FormCalc, and that Cuba's original definition is slightly more involved.) 

• A partition viewer has been added. It allows to view the partitioning of the integration region performed by 
Suave, Divonne, and Cuhre (Vegas does not tessellate the integration region). To use the partition viewer, set 
the verbosity level to 3 and pipe the output of your program through the partview program, as in 

myprogram I partview 12 13 
Each pair of numbers on the command-line refers to a hyperplane to display, i.e. in the above example the 1-2 
and 1-3 hyperplanes are shown. See Fig. fl] for a screenshot. 
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Figure 1: A screenshot of the Cuba Partition Viewer. 



Cuba possesses also a Mathematica interface. The MathLink executables are loaded e.g. with Install ["Vegas"] 
and make the Cuba routines available in a way almost like NIntegrate. 

Mathematica C 
MathLink 



Vegas [/, 

, L_ 



.] 



Compile 



T 

integrand / 



{xi,x 2 ,-.-} 
{fi,h,-..} 



void Vegas (. . . ) 



request samples 



The integrand stays in Mathematica at all times (see figure above), which means that one can integrate functions 
which are not simple to implement in C or Fortran, e.g. Cuhre [Zeta [x y] , {x , . 2 , . 3} , {y , . 4 , . 5}] . 



4. Public code-generation functions 

FormCalc's code-generation functions, hitherto used only internally, have been reorganized for greater user- 
friendliness and are now publicly available. They can be used to write out an arbitrary Mathematica expression 
as optimized Fortran code. The basic procedure is very simple: 

1. handle = OpenFortran I" file. F"^ 
opens file.F as a Fortran file for writing, 

2. WriteExpr [handle, {var -> expr, . . . }] 

writes out Fortran code to calculate expr and store the result in var, 

3. Close [handle'] 
closes the file again. 



0604 



The code generation is fairly sophisticated and goes well beyond merely applying Mathematica's FortranForm. The 
generated code is optimized, e.g. common subexpressions are pulled out and computed in temporary variables. 
Expressions too large for Fortran are split into parts, as in 

var = parti 

var = var + part 2 

If the expression is too large even to be sensibly evaluated in one file, the FileSplit function can distribute it on 
several files and optionally write out a master subroutine which calls the individual parts. 

To further automate the code generation, such that the resulting code needs few or no changes by hand, many 
ancillary functions are available, e.g. CommonDecl writes out common-block declarations for a given list of variables. 

5. FeynHiggs interface and NMFV 

FormCalc's initialization code for the MSSM, mssm_ini .F, now contains an interface to FeynHiggs 0, for 
the computation of the Higgs masses. Since in particular the light Higgs-boson mass receives significant radiative 
corrections, its precise determination is phenomenologically important. From the user perspective, a preprocessor 
flag governs the choice of Higgs masses: 

• #def ine HIGGS_MASSES TREE 

use the tree-level Higgs masses, 

• #def ine HIGGS_MASSES SIMPLE 

use a simple one-loop formula, 

• #def ine HIGGS_MASSES FEYNHIGGS 

invoke FeynHiggs 2.2 |8| - this is the most precise determination, 

• HIGGS_MASSES undefined 

use the FeynHiggsFast approximation Q - quite precise, but valid only for real SUSY parameters. 

Also for the MSSM, non-minimal flavour violation (NMFV) can be enabled. In the usual, minimal, setup, there is 
no flavour-violation beyond the Standard Model's CKM effects, i.e. left -right mixing is independent for each squark 
flavour. In contrast, NMFV means full 6x6 mixing among both up-type squarks (u\, U2, c%, £2, £1, £2) and down-type 
squarks d 2 , h, h, h, h), see e.g. 0. 

On the technical side, the diagrams have to be generated using the FVMSSM.mod model file in FeynArts, and in 
the numerical evaluation the preprocessor flag FLAVDUR_VIDLATIDN has to be defined. Once this flag is enabled, the 
6x6 mass matrix is accessible as LambdaSf and in addition to the usual squark masses and mixing matrices, MSf 
and USf , there exist MASf and UASf containing the corresponding NMFV quantities. 

6. Shell scripts 

FormCalc 4.1 includes a few useful shell scripts: 

• sf x packs all source files in the directory it is invoked in into a mail-safe self-extracting archive. The archived 
code in particular does not need FormCalc to compile or run. 

• turnof f switches off (and on) the evaluation of certain parts of the amplitude, which is a handy thing for 
testing. For example, "turnof f box" turns off the boxes (actually, all parts of the amplitude with 'box' in 
their name), turnof f without arguments then restores all modules. 

• pnuglot produces a high-quality plot in Encapsulated PostScript format from a data file in just one line. 
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• submit automatically distributes a parameter scan on a cluster. The available machines first have to be declared 
in a file .submitrc, e.g. 

# Optional: 'nice' to start jobs with 
nice 10 

# Pentium 4 
pcl301 
pcl305 

# Dual Xeon 
pcl247b 2 
pcl319a 2 

After that the command line, typically something like "run uuuu 500,1000", simply has to be prefixed by 
submit, i.e. "submit run uuuu 500,1000". 



7. The HadCalc frontend 

HadCalc is a new frontend for FormCalc, i.e. it uses the generated Fortran code with a custom set of driver 
programs. It was written by Michael Rauch who also maintains it independently from FormCalc. 

HadCalc automates the calculation of hadronic cross-sections. It automatically performs the convolution with 
PDFs and allows various cuts to be applied. Cross-sections can be computed either fully integrated or differential 
in invariant mass, rapidity, or transverse momentum. HadCalc operates either in batch mode, like FormCalc, or 
interactively, which allows the user e.g. to play with parameters. 

The program is not (yet) public. It can currently be obtained on request from (mrauch@mppmu.mpg.de). 



8. Summary 

FormCalc 4.1 is the current release of the FormCalc package. It has many new features enhancing performance and 
user- friendliness. The new public Fortran code-generation functions should be interesting also to people who do not 
use FormCalc to compute Feynman amplitudes. The package is available as open source and stands under the GNU 
Lesser General Public License (LGPL). It can be obtained from the Web site http://www.feynarts.de/formcalc. 
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